Oscillations of the superconducting critical current in Nb-Cu-Ni-Cu-Nb junctions 
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We report on experimental studies of superconductor-ferromagnet layered structures. Strong 
oscillations of the critical supercurrent were observed with the thickness variation of the ferromagnet. 
Using known microscopic parameters of Ni, we found reasonable agreement between the period of 
oscillations and the decay of the measured critical current, and theoretical calculations. 



The interplay between superconductivity and ferro- 
magnetism is an oMj^ubject which was studied exten- 
sively over decadesErEl. The most striking effect in such 
systems is the formation of the so called tt phase junction 
in a superconductor-ferromagnet-superconductor (SFS) 
structured. One of the manifestations of the tt phasG||i6 
a non-monotonic variation of the critical temperatureo'Q, 
Tc, with the variation of the ferromagnetic layer thick- 
ness, d. Oscillatory dependence of Tc vs. ^(L of SF 
multilayers has begii— observed by few groupsQEl, how- 
ever, other groupsBllj reported monotonic behavior in 
similar structures. Interpretations of these experiments 
suggested the tt juncti(m_raechanism as well as soroe 
other origins of the effcctEllljO. The recent observationEj 
of non-monotonic behavior of the critical current as a 
function of temperature in weak-ferromagnetic layer of 
Cua;Nii_a; between two Nb layers is considered as an 
unambiguous proof of the tt phasfj-fajytaation. Another 
interesting theoretical predictionl3fiEJli3 concerns non- 
monotonic behavior of the critical current as a function of 
the thickness of the ferromagnetic layer d. According to 
the above predictions, the critical current Ic is expected 
to oscillate and decay as d is increased. To the best of 
our knowledge, such a behavior has not been reported so 
far. 

In this paper, we present the experimental evidence of 
oscillatory behavior of the critical current vs. thickness 
variation of ferromagnetic Ni layer. We also show a 
sonable agreement between our data and the theoriea3'E£l 
in the appropriate limit of Ef.x » ^/t >> kBTc- Here, 
r and E^x are the electron relaxation time and the ex- 
change energy of the ferromagnet, and Tc is the critical 
temperature of the superconductor. 

We have studied temperature and thickness depen- 
dence of the critical current in Nb-Cu-Ni-Cu-Nb junc- 
tions. Additionally, we have studied the thickness de- 
pendence of critical current of similar junctions without 
a Ni layer. The junctions with lOxlO/im^ area were fab- 
ricated with the standard photolithography technique. 
The process contained three stages of lithography: liftoff 
of the bottom Nb-Cu layer, liftoff of the variable thick- 
ness Ni or Cu layers, and liftoff of the top Cu-Nb layer. 
Nb films were sputtered using a magnetron gun and in 
situ covered with the Cu layer by thermal evaporation, 
for preventing the Nb oxidation. The ferromagnet lay- 
ers of Ni were e-gun evaporated in a separate vacuum 
chamber, and subsequently covered in situ by Cu. It is 
important to emphasize that all samples were prepared 



simultaneously. The variation of Ni thickness for Nb-Cu- 
Ni-Cu-Nb junction and variation of Cu for Nb-Cu-Nb 
junctions was obtained by a specially designed shutter, 
which exposed the samples in sequence, so that every 
sample was exposed to the evaporating Ni or Cu for ad- 
ditional fragments of time. This method guaranteed that 
all the interfaces between each layer in our multilayer 
structure are identical, and the only difference between 
the samples is their Ni or Cu thickness. The thickness of 
each Nb layer was 2000A. The total thickness of the Cu 
in the Nb-Cu-Ni-Cu-Nb junctions was 2400A and the Ni 
thickness varied from lOA to 90A. 

In the other set of junctions, namely Nb-Cu-Nb, the 
Cu thickness varied from 5000A to 10000 A. The struc- 
ture of these junctions is shown schematically in the inset 
ofFigj|. 

Typical resistance of both types of junctions was 
around 100/xri, independent of both Ni and Cu thick- 
nesses. The critical temperature of the Nb was about 
8.5K. The measured resistivity |£if Ni layer p = Slfiflcm^ 
together with the known valuetj pi = 1.5 • 10~^^flcm^ , 
allowed us to establish the mean free path / = 48A. 
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FIG. 1. dV/dl of one of the samples. The critical current 
was defined as the onset of the peak. In the inset; schematic 
picture of a Nb-Cu-Ni-Cu-Nb junction. 

The measurements were performed in "^He cryostat in 
the range from 6.5K down to 1.5K. The critical current 
was measured by passing a DC current with a small AC 
modulation through the sample. The AC voltage which 
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appeared above the critical DC current was picked up by 
lock-in amplifier operated in a transformer mode. 

A typical differential resistance of both types of 
junctions vs. DC current is plotted in Fig.|l|, which 
clearly indicates on the sharp onset of dissipation at 
I = Ic- First, we show results of the Nb-Cu-Nb set 
(Fig.^). As expectedllll, the critical current Ic normal- 
ized to the length of the junction L, decreases expo- 
nentially with L. From the decay of the critical cur- 
re nt we found the thermal length of the the Cu, Lt — 
y^hD/2'KkBT = I6OOA at 4.2K, and the diffusion con- 
stant D — 860c?7i^/s. 
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FIG. 2. Tlie normalized critical current Ic/L of a 
Nb-Cu-Nb junctions as a function of the Cu thickness L. The 
straight line represents the theoretical decay. 

Very different thickness dependence was found in the 
Nb-Cu-Ni-Cu-Nb junctions. Fig.|| shows the thickness 
dependence of the critical current in these junctions at 
T = A.2K. In spite of the large error bars, the non- 
monotonic variation of the critical current is quite evi- 
dent, namely, the deviations of the data from the expo- 
nential decay surmounts by far the uncertainty of each 
measured point. 

A further, and even a stronger evidence for the oscil- 
latory behavior is provided by temperature dependence 
of the critical current. Fig.^ shows a family of Ic vs. 
T curves for all Ni thicknesses, which are normalized to 
their values at A.2K . We define the slope ax of temper- 
ature variation, namely ar = d{lnIc(T))/dT, and plot 
these values as a function of d in Figjl (squares). Oscil- 
lations of ax are very prominent, and are in anti-phase 
with the oscillations of the critical current. Unlike the 
critical current, which had quite large experimental error 
bars, the slope of the temperature dependence had error 
bars of only few percents. Such a behavior of ax as a 
function of d is intimately related to the variation of the 
critical current oscillations amplitude with temperature. 
Although the amplitude of the oscillations decreases at 
high temperatures, their relative value increases, and on 



a semi-logarithmic plot of Fig.||, the oscillatory behavior 
of Ic is more pronounced at higher temperatures. 
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FIG. 3. Critical current of the Nb-Cu-Ni-Cu-Nb junctions 
as a function of the Ni layer's thickness d at 4.2K (cir- 
cles). The dependence of the slope ar on d is represented 
by squares. Both the dashed and dotted lines are only for 
guiding the eye. 

We would like to start our discussion with what is well 
known for superconductor-normal-superconductor (SNS) 
junctions. The dependence of critical current of such 
junctions on the thickness of N layer is exponential, 
Ic a e~^/^^ for L > LtiH. In the opposite limit, where 
L < Lt, Ic approaches a constant value. This behav- 
ior has been obs er . ve d in few experiments for different 
material systemsli3iij, including our sample presented in 
Fig.|. 
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FIG. 4. Critical current as a function of temperature of the 
Nb-Cu-Ni-Cu-Nb junctions for different thicknesses of the Ni 
layer. 

Non-monotonic behavior of the critical current is not 
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expected theoretically, and had never been observed ex- 
perimentally in SNS structures. Since the thickness of 
Ni is the only parameter which is varied in our Nb-Cu- 
Ni-Cu-Nb samples, it is reasonable to assume that the 
oscillations observed in the data are due to the presence 
of Ni in the structure. 

As mentioned in the introduction, oscillatory behavior 
of the critical current vs. the thickness-oif the ferromag- 
netic layer is predicted theoreticallyBa^tj. The origin of 
these oscillations is the phase shift acquired by electron- 
hole Andreev particles upon entrance into the ferromag- 
net, due to their different spin orientations. Several ex- 
pressions have been derived for the critical current in var- 
ious limits of the strength of Eex , thickness of the ferro- 
magnet d and disorder. Since in our experiment we have 
determined only the magnitude of the critical current /c, 
the formulae below will be written for the absolute value 
\Ic\- For the clean, thin and strong ferromagnetic layers 
I > hv f / Eeand, the Critical current should vary with the 
thickness asEl: 
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FIG. 5. Critical current as a function of the thickness of 
the Ni layer at different temperatures. 
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The dependence of absolute value of the critical current 
is shown on Fig.^ by dashed line. In the opposite limit, 
namely / < d, the following dependence was derivedEI 
(dotted hue in Fig.^: 



/e ~ ye-y |sin(y + 7r/4) | y = yJhD/^E,, 



(2) 



Another expression for the critical current in the limit 
I > fiVf/Ef,^, is given bytj: 



ReY: 
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(3) 



The curve given by Eq. || in Fig.^ follows the decay 
and the phase of the oscillations of the measured Ic closer 
than the other two curves for large d. Furthermore, we 
estimate hvf/Eex — 17 A which is smaller than the mea- 
sured mean free path I = 48A, and therefore, the validity 
of using Eq. |^, / > hvf/Eex, is justified. However, since 
Eq. H is valid only for d > I, the data points should fol- 
low Eq. H for d < 40A. We thus conclude that our data 
is consistent with Eq. ^ and Eq. ^ in their appropriate 
limits. Therefor, we give the fit of Eq. and Eq. | to 
our data for four temperatures in Fig. |^. 



where u>n = TrTkBi^n -I- 1) is the Matsubara frequency, 
n is an integer number, k^^ ^ {I + 2|cij„|r/?i) — 2iEg^xT/h, 
fj, = cos9, 9 is the angle between the momentum and the 
normal to the SF interface, and A is the order parameter 
in the superconductor. The solid line of Fig. ^represents 
the above expression. In fitting cuir data to the expres- 
sions for all limits, we have usedE3 Vf = 2.8xl0^m/sec 
and / — 48A (based on the resistivity measurements of 
Ni). Therefore, the only fitting parameter apart from 
the numerical prefactor was the strength of the exchange 
interaction Eex- 

The periodicity of oscillations in Eq. ^ Lose 
nhvf/Eex — 54A fits the best our data when Eex = 
107 ± 3meV. Note, that Lose is twice larger than the 
periodicity observed in Fig.^, due to the absolute value 
taken for le- In order to have similar periodicity in Eq. 
^, we have used weaker exchange field Eex = SGmeV^ 
All values are close to the recently reported valueE2l 
Eex = 115meV. 



< 

E 




d[A] 



100 



FIG. 6. Critical current of the Nb-Cu-Ni-Cu-Nb junctions 
at 4.2K, and the theoretical fits. The dashed line represents 
Eq. 1, the dotted line Eq. 2, and the solid line Eq. 3. 
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FIG. 7. Critical current of the Nb-Cu-Ni-Cu-Nb junctions 
at 1.5K, 4.2K, 5.2K and 6.5K, and the theoretical fits accord- 
ing to Eq. 3 (solid line) and Eq.I (dashed line). 

The oscillations of ut (Fig. which slightly smear 
the oscillations of Ic (Fig. |^), are not accounted by any 
of the mentioned theories and have only a minor effect 
on the overall quality of the fit presented in Fig. |^. A 
possible origin of this effect could be related to the pres- 
ence of normal metal on both sides of the ferromagnetic 
layer. The structure we have used contains layers of Cu of 
a combined thickness L ~ 2400A, which is much larger 
than the thickness of the Ni layer. It is obvious that 
at high temperatures {Lt < L) the contribution of An- 
dreev particles propagating via short diffusion path is 
predominant, since paths which are much longer than L 
do not contribute to the supercurrent. At low temper- 
atures, where Lt > L, the contribution of longer paths 
is not negligible. The latter implies that the impact of 
the trajectories crossing Ni layer several times is higher 
at low temperatures. The multiple crossing of the ferro- 
magnetic layer could produce additional harmonics of the 
basic frequency of the oscillations, which will smear the 
non-monotonic behavior. This multiple crossing has not 
been considered in the existing theories for S-F-S junc- 
tions, and therefore might be responsible for the observed 
temperature dependence of the oscillations in S-N-F-N-S 
junctions (Fig. |). 

In summary, we have observed the oscillations and the 



decay of the critical current in Nb-Cu-Ni-Cu-Nb junc- 
tions upon the increase of the thickness of the Ni layer. 
We found a reasonable agreement with the recent theo- 
retical calculations in the appropriate limit. The puzzling 
temperature dependence of the oscillations' amplitude re- 
mains unexplained. 
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